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o : 

• The conformational behavior of a polymer in a critical binary solvent confined in a porous medium 

OA ' is studied. The size of the polymer in bulk, which is mainly governed by the correlation length of the 
solvent density fluctuations, depends on the proximity to the critical point of the binary mixture. 

^ ' We find that in contrast to the bulk behavior, the conformational properties of the polymer in a 

I porous medium depends strongly on the pore size. The latter controls the correlation length of the 

. solvent density fluctuations and thus determines the polymer size. 

cJi ■ I. INTRODUCTION AND MODEL 

o : 

^ ; ' The conformational behavior of polymers in a critical (binary) solvent is mainly governed by two length scales - the 
d chain length and the correlation length of the solvent density fluctuations [1,2]. In particular, the size of a polymer 
' in bulk is determined by the proximity to the critical point of the solvent mixture. Far away from the critical point, 
I I the size of a polymer in a binary mixture of good solvents simply scales as i? ~ iV"^/^, where N is the number of 
' O monomers. Close to the critical point, however, when the correlation length of the solvent density fluctuations is 
^ comparable to the chain size, the polymer collapses to form a globule of size R ^ N^/^; a collapse transition close 
to the critical point is believed to be driven by the solvent density fluctuations which induce attractive interactions 
between different parts of the chain. At the critical point itself, when the solvent density fluctuations act on a much 
larger scale compared to the polymer size, the chain regains its original size, that is, R ~ N'^^^. These effects first 
studied by Brochard and De Gennes [2] using simple scaling arguments were later confirmed by simulations of Magda 
et al. [3], field theoretical and variational studies [1,4,5] and experiments [6-8]. 

The fact that the solvent density fluctuations are instrumental in determining the bulk behavior of polymers in a 
critical binary solvent raises an important question when such a polymer mixture is present in a restricted geometry, 
^\ • say a porous medium. Since the pore size affects the correlation length of the solvent density fluctuations and thus 
' the polymer size, the question as to how a given confinement controls the conformational behavior of the polymer 
\l is an important one. This paper addresses this issue, which, we believe, has its implications in certain experimental 
situations like chromatography. 

We begin by presenting a brief review of the bulk behavior — a self-avoiding polymer chain immersed in a binary 
solvent mixture, and describe it using the Edwards Hamiltonian formalism. The idea is to determine the effective 
potential exerted on a chain under various conditions of temperature and fluid composition. The Hamiltonian that 
. describes the system of interest is given by 
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where (3 = l/k^T, ks is the Boltzmann constant and T the absolute temperature; R(s) is the chain variable, s the 
curvilinear coordinates along the chain (contour variable) ; the chain consists of N segments each of size b. The first 
term accounts for the connectivity of the chain and represents the chain entropy due to its elasticity. The second term 
results from two body interaction between monomers, where v is the strength of the bare monomer-monomer excluded 
volume interaction, which is assumed to be short ranged. The third term represents the interaction between the chain 
and the fluid in which Va- is the strength of the short range excluded volume interaction between the monomers and 
the different species of the fluid, vi ^ V2, i.e. the solvent quality is slightly different. In essence, it means that the 
chain is likely to be surrounded by solvent 2, when vi > V2- /37if describes the Hamiltonian of the fluid, the form 
of which needs to be specified. At this point we introduce a collective density field (order parameter) for the fluid, 
which is given by 
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c(r)=^<5(r-rO, (2) 



which allows to rewrite the Hamiltonian in the simple form 



r-N 

+ {vi-V2)y2 rfsexp[-ikR(s)]c(k)+/3W/(c(r)), (3) 
k -^0 

where PTi.{{c) can now be viewed as 

pHMr)) = I d'r Q(|Vc(r)|2 +rc2(r)) + A^4(r)^ . (4) 

The parameter t describes the distance from the critical point and A is a coupling constant. Far away from the critical 
point the quartic term in Eq. (4) can be neglected, which amounts to A = 0. The latter is a Gaussian approximation, 
which allows us to integrate out the solvent degrees of freedom to produce an effective interaction potential between 
different parts of the chain [1]. The effective Hamiltonian is given by 

mes =^j^ {^^) "^'^^C C '^''''^^^ exp[-ik(R(s) - R(s'))], (5) 

where 

In writing the above equation, we have defined A?; = vi — V2 ■ Since the (mean field) correlation length of the solvent 
density fluctuations is defined by ~ b/y/r, the form of Eqs. (5) and (6) suggests that in certain temperature 
regimes depending on the strength of interaction between the chain and the fluid, the effective potential can become 
attractive (negative) to induce a chain collapse. In other words, depending on the correlation length of the solvent 
density fluctuations and At;, the effective solvent quality can become poor. Therefore, before we discuss the problem 
of a chain in a restricted geometry, it is useful to briefly review the behavior of a chain in a nanopore under poor 
solvent conditions. 



II. NANOPORES 



A. Chains in Nanopores 



Let us first review some of the scaling results of chains in nanopores. The diameter D of the tubes is assumed to 
be smaller than the natural chain radius Rg ~ A'''', where in three dimensions u = 1/2 for theta-solvents and i/ w 3/5 
for good solvents. 

We consider flexible polymer chains confined in pores of cylindrical shape with neutral walls such that there is 
no adsorption between the polymer and the pores; the volume available to the solvent and the solute is limited by 
well defined boundaries. For the case of a single chain in a capillary of diameter D, there are two distinct behaviors 
depending on the dimensionless ratio Rp/D, where Rp is the Flory radius of a polymer in bulk. If Rp < D, we 
have a conventional bulk solution. If Rp > D the chain is deformed and is confined in a cylindrical tube of diameter 
D <C Rp- Since D ^ h, the chain still retains some lateral wiggling. We assume that the tube wall repel the chain 
strongly so that there is no trend towards adsorption. The length of the tube occupied by the chain is R\^. The 
behavior of the chain can mainly be understood with the blob picture. The size of the blob is solely determined by 
the tube diameter. The parallel size of the chain is then made out of rib blobs of diameter D such that R\\ — ni,D. 
The only information about the quality of the solvent is then in the number of blobs themselves, which is given by 
rib = N/g, where g is the number of monomers inside a blob. The latter is simply given hy g = {D/hY^" , and the size 
of the chain along the tube is ~ bN{b/ D)^^~"^/'^ . Note that the same results can be derived by using a free energy 
as I3T = Rj/Nb'^ + pU, where pU ~ b^N'^/D'^R\\ for good solvents and I3U ~ b^ / {D'^ R\\f for theta-solvents. 
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As discussed in the last section, the solvent density fluctuations can induce effective attractive interaction between 
different parts of the chain segments leading to a chain collapse. Since this amounts to the effective solvent quality 
becoming poor, it is worthwhile to consider the case of a polymer in a poor solvent confined in a nanopore. The free 
energy of such a system is given by 

Rf. iV2 ^3 

where v ^ —tit' and w ^ are the two and three body interactions respectively. Here the reduced temperature t 
describes the distance from the theta temperature, i.e., t = \T — 0|/0. Let us first discuss the influence of the tube 
geometry on the O - point. To estimate the 9— temperature, we employ the usual Ginzburg criterion. The 6— point 
for a flnite chain length in a restricted geometry cannot be defined by simply setting the second virial coefficient 
equal to zero. Instead, a consistent expression can be determined by the balance between the chain entropy and the 
attractive second virial term. In particular, when t is sufficiently large (the solvent quality is poor enough), the first 
term in Eq. (7), which represents the entropic elasticity, is smaller than the attractive interaction due to the second 
term. Thus the Q— point for chains inside the pore can be determined by comparison between the entropic term and 
the negative second virial term. This yields the following condition for the G- point: 

t:^{D/bf{R\i/bNf, (8) 

where the estimate is sensitive to the unperturbed chain size . The usual Ginzburg criterion uses the ideal chain 
size R oc bN^/^, but this yields the unphysical result tc — {l/y/N){D /bN^/'^Y since it suggests an increase in t with 
the increase in the pore diameter. Moreover, the chain length cannot play a significant role since the collapse in the 
pore takes place only on the length scales less than the pore diameter, i.e., the relevant blob size. Therefore, it is 
physical to use the Gaussian value for R\\ inside the pore, which is simply given by i? ~ bN{b/D), as suggested by a 
simple blob argument. As a result, the value for the shift of the G - temperature is given by 

* = i (9) 

The above equation shows significant dependence of the G-temperature on the pore size. It can be rewritten as 

9 = 6'buik + const (10) 

which is only valid for D < b^N, i.e., when the pore diameter is smaller compared to the Gaussian chain size. When 
D becomes comparable to the chain size, then the shift of the G - temperature is determined by t ~ b/RGauss, which 
is simply t oc I/a/ZV, as it must be the case for the collapse of polymer chains in bulk. The geometry dependence 
of the G-tcmpcrature has recently been studied for a polymer confined between two parallel plates (slit geometry) 
by using numerical methods [12]. For completeness, we would like to remark that naive scaling for a simple slit 
geometry predicts tc,siit = (b/D)^, which suggests a logarithmic dependence, tc.siit oc \og{D/b), but this needs further 
investigation. 

For t > tc the chain becomes significantly contracted from its theta size. The stretching term can then be neglected, 
and the balance between the attractive and repulsive third body potential yields the chain size 

b./b^' 



i^„-iV^-J. (11) 

The density inside the tube is simply given by p ~ ~ t/b^, which corresponds to the density of a unconstrained 

globule and also defines the thermal blobsize = b/t. On length scales smaller than the chain statistics are 

1/2 

unperturbed by the volume interactions and are that of a random walk, i.e., ~ bg^i . The parallel extension of the 
chain can, therefore, be expressed by comparing the blobsize and the diameter D, 

R\\ - bN^-^ (i) ■ (12) 



D 

Using the above expression, the size of the polymer confined in a tube of diameter D can be summarized for three 

distinct cases: When D < £^t, the chain experiences theta conditions; when D > ^7-, the chain inside the tube finds 
itself in a poor solvent condition. It collapses to form a globule of size governed by t. The chain starts to shrink when 
D = ^T, starting from the size Ru ~ bNt, i.e., a linear arrangement of the thermal blobs. 
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B. Critical fluids in Nanopores 



The problem of phase transitions in restricted geometry has been considered mostly in the case of (two dimensional) 

wetting. In a nanoporc the situation is very different. The problem of a spherical model in a sphcrocylindcr has been 
studied by Cardy [13], who shows using conformal mapping that the correlation length at the critical point can be 
written as 

= (13) 

where A = d — 2 + 7] is a universal amplitude; it depends only on the space dimension d and the critical exponent 
T]. this results implies also that the amplitude and so ^ does not depend directly on the value of the fluid coupling 
constant A, which was introduced in eq. (4). Physically, this implies that the criticality plays a significant role only 
on the scale of the tube diameter D and fluctuations always stay finite. Therefore, the effects of the solvent density 
fluctuations on the conformational behavior of a polymer chain, which is dissolved in a bicomponent solvent present 
in a porous medium, are mainly local. These effects are studied in the next section. 



C. Polymer chain and fluids in nanopores 



As discussed in the previous section, when the chain dissolved in a bicomponent solvent is confined in a nanoporc, 
the problem is effectively one dimensional. To see this we can rewrite all spatial vectors in their Cartesian coordinates, 
that is 

fj,=x,y^z ^ ^ ^ u u f_L—x,y,z 

I'M 

+ {vi-V2)Y,j dseM-i E A:^fi/.(s)]c(k)+/3H/(c(r)). (14) 

Moreover, when the constraints due to the cylindrical geometry of the tube are imposed the effective potential 
remains structurally similar. Because of the confinement, the vector k can be splitted into its Cartesian components, 
i.e., k = {k\[,k_\_,k±_), where fcy corresponds to the length scale along the (cylindrical) tube. The component k±_ is 
limited by the tube wall and can be replaced by its upper limit 1/D. This naive estimate suggests that the effective 
potential is given by 

v{^.) = v (15) 

Within this mean field model, the new correlation length can the be identified as 

b%l. = r+{^^)\ (16) 

It is interesting to note that the scaling of the correlation length agrees with the exact result of Cardy [13] at r = 0. 

To estimate the size of the chain within the classical scaling, it is sufficient to discuss the limit = 0. The effective 
potential can then be casted into a Flory form, which yields together with the Gaussian stretching term a reasonable 
estimate for the size of the chain. The effective potential is given by 

v = v- ^ ' „ ■ (17) 

For T = 0, i.e., at the critical point this equation is similar to the corresponding mean- field equation in the bulk, Eq. 
(6) in the low wave vccotr limit. For geometrically unrestricted systems D oo, the correlation length grows then 
with the pore diameter, because here we have oc D. At the critical point, i.e., r = 0, the effective monomer-monomer 
potential v becomes negative when the tube diameter is sufficiently large: 



(Wl - V2) 
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Since this amounts to the effective interaction between the chain segments being attractive, the chain is expected to 
collapse inside the tube. In contrast to the bulk behavior, where in a certain temperature regime (close to the critical 
point) the chain is always collapsed, the collapse in a confined geometry depends purely on material properties — the 
strength of the polymer - fluid interaction, the chain excluded volume, and the diameter of the tube. The geometry, 
therefore, plays an important role in the collapse of a polymer chain both in a critical fluid and a fluid with a flnite 
correlation length. 



D. The size of the chain 



The Hamiltonian as given by Eq. (5) can be transformed into a Flory free energy by simple dimensional analysis 
[11]. It has to be modified such that the 3-dimcnsional Dirac function is anisotropic and the lateral dimensions are 
determined by the diameter D of the nanopore. Thus, we can estimate the excluded volume to be (5(R(s) — R(s')) oc 
where is the chain size parallel to the pore. The interaction parameter between the monomers is 
determined by the value for the effective chain potential ^(fcii = 0): 

/3^=^+^(^ll=0)^ (19) 
Minimization of the above free energy yields the size of the chain in the nanopore: 



D 



The above expression is valid as long as the effective potential is positive. At the critical point r = 0, we can use 
Cardy result for d = 3 and to get 

where we have rescaled the interaction parameters by the size of the Kuhn length v ^ It' and introduced = (''^i^W2)^- 
This results holds as long as the effective potential is positive, which means that the pore size is small enough. For 
larger pore sizes as 

D > 7 (22) 

the effective potential becomes negative; the critical fluid then resembles a poor solvent. The size of the chain can 
then be estimated by taking into account the three body interaction in the tube. It lateral extension is then given by 

1 /b-^ 



Therefore, the chain is always in a stretched conformation, as has already been seen from the scaling analysis of the 
behavior of chains in solvents of different quality. 



III. APPLICATIONS AND DISCUSSION 



The behavior of chains immersed in a critical fluid confined in a porous media is substantially different from the 
bulk case. As mentioned in the Introduction, a chain in bulk undergoes a conformational transition, i.e., N^^^ — > 
JV^/^ N^/^^ on approach to the critical point of the fluid. In a restricted geometry, however, the correlation length 
of the fluid stays finite and is mainly determined by the typical size of the tube. The latter implies that the pore 
diameter can select chains according to their chain length — for a given pore diameter D, short chains (of size much 
less than D) arc not restricted by the confinement; they can midcrgo a collapse transition on approach to the critical 
point, as is the case in the bulk. Sufficiently long chains (of size much larger than D), on the contrary, remain in a 
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stretched conformation. As a result, the variation in the fluid temperature, that is, approach towards its "criticality" 
can select chain conformations according to their polymerization degree. 

To see this let us consider a dilute solution of chemically identical chains having degrees of polymerization iVi and 
N2 dissolved in a good solvent such that Ni » N2. In terms of the present model, a good solvent quality implies 
T > — {b/D}^ (i-G- the temperature not too close to the critical point where r = 0). When the pore size is such 
that D N2^^, only long chains of size N^^^ will be confined. The confinement free energy for the long chain is given 

by 



5/3 



f3F ^ N, ^-j , (24) 

which is estimated as the number of blobs times fc^T. The free energy for the non-confined chain is roughly given by 

Analog 7V2 ~ (7- 1) log ^^2. (25) 
where 7 is the partition function critical exponent. The comparison of the two free energies yields 



D /A^i^^/^ 



b \N2 



(26) 



which provides a simple scaling for the separation of chains according to their chain length. It is to be noted that in 
writing the above expression we have ignored the logarithmic dependence in Eq. (25). 

Eq.(9) suggests another interesting possibility — the chain length dependence of the second virial coefficient allows 
us to determine the temperature away from the critical point where the separation between the collapsed short chains 
and the more stretched long chains occur: 

r* = TT477^-(i^X- (27) 



l + {b/D) \D 

For chains whose size is less than that of the pore diameter, the above condition suggests 

= (28) 

Thus by tuning the distance from the critical point, that is r, the chains may adopt different (globular) conformation 
depending on their chain lengths. This provides a means for separating chains by changing the poor solvent quality 
(through the change in the criticality of the fluid) in porous media. 
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